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Abstract. New physics near the TeV scale could modify neutrino-matter interactions or 
generate a relatively large neutrino magnetic (transition) moment. Both types of effects have 
been discussed since the 1970's as alternatives to mass-induced neutrino flavor oscillations. 
Nowadays, the availability of high-statistics data makes it possible to turn the idea around and 
ask: How well do the simple mass-induced oscillations describe solar neutrinos? At what level 
are the above-mentioned nonstandard effects excluded? Can we use solar neutrinos to constrain 
physics beyond the Standard Model? These notes review the sensitivity of the present-day 
solar neutrino experiments to the nonstandard neutrino interactions and transition moment 
and outline progress that may be expected in the near future. Based on a talk given at the 
Neutrino 2006 conference pQ. 
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1. Standard LMA solution: basic features 

The most basic experimental fact about the neutrinos from the Sun is that the electron neutrino 
survival probability, P^f" = P(y e — > v e ), is measured to vary as a function of the neutrino 
energy. At the high end of the spectrum (E v > 6 — 7 MeV) the SNO [2] and Super-Kamiokande 
[3] experiments have established that P^ d is about ~ 34 ± 3%. The gallium experiments [3], 
however, which are sensitive to both high- and low-energy neutrinos, see a higher survival 
probability: the measured rate is 74 ± 7 SNU, whereas the standard solar model prediction 
[5] (before oscillations) is 1311Jq. 

This simple fact has highly non-trivial implications. Indeed, this behavior is not "generic" 
for mass-induced oscillations, even when they combine with the MSW [U \7\ matter effect. A 
priori, one might have expected solar neutrinos to be in one of these regimes: 

• matter dominates at the production point, on the way out of the Sun neutrino flavor evolves 
adiabatically — > constant suppression (regime 1); 

• matter dominates at the production point, on the way out of the Sun neutrino flavor evolves 
non-adiabatically — > vacuum oscillations 

— vacuum oscillation length <C 1 a.u. (astronomical unit) — > oscillations average out — > 
constant observed suppression (regime 2); 

— vacuum oscillation length S> 1 a.u. (astronomical unit) — * no time to oscillate — * no 
suppression (regime 3); 

• vacuum oscillations dominate everywhere, matter effects negligible even in the center of the 
Sun — ► oscillations average out — * constant observed suppression (regime 4). 



The observed energy-dependent P^ d then implies that solar neutrinos are in one of the several 
"special regimes": the transition between regimes 1 and 4 (Large Mixing Angle - LMA - 
solution); the transition between regimes 2 and 3 (vacuum/quasi- vacuum [8] oscillation solution); 
the transition between regimes 1 and 2 (Small Mixing Angle - SMA - solution); the regime where 
the density in the Earth is close to resonant, so that the flavor regeneration in the Earth is large 
(the LOW solution). These solutions were known for many years, in particular all four were 
allowed as recently as 2000, see, e.g., [9]. 

We now of course know that only the LMA solution survives. Let us consider the survival 
probability P^e under the (a posteriori justified) assumption that the oscillations take place 
between just two eigenstates. One easily obtains that during the day time 



p£ d > 2 » = cos 2 # cos 2 9 + sin 2 6 Q sin 2 i 



(1) 



The probability of finding the neutrino in eigenstate 1(2) is cos 2 #0 (sin 2 9®), where 9® is the 
mixing angle at the production point; in turn, the probability of detecting the neutrino already 
in eigenstate 1(2) as v e is cos 2 (9(sin 2 9). The key physical ideas here are that the evolution 
is adiabatic (no level jumping) and incoherent (interferences between 1 and 2 disappear upon 
integration over energies for Am? > 10 -9 — 10 -8 eV 2 \W\ [8] and over the production region). 



The angle 9® is determined from the oscillation Hamiltonian H to t = H V£LC + H m£tt , where 
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Here A = Am 2 / (AE V ) and Am 2 is the mass splitting between the first and second neutrino 
mass states: Am 2 = m| — m\. The two limiting values are 9® = 9 (H tot is dominated by H vajC ) 
and 9® = ir/2 (H to t is dominated by H ma , t )- The probability P^f then varies from cos 4 6 l + sin 4 9 
(= 1 — (1/2) sin 2 29, averaged vacuum oscillations) to sin 2 9. 
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Figure 1. The v e survival probability and day/night asymmetry for the LMA solution. 



The transition from one regime to another occurs when \[2Gpn e ~ 2 A at the production 
point. To accommodate the data on P ee , this transition must occur right in the middle of the 
solar neutrino spectrum, implying Am 2 ~ a few x 10 -5 eV 2 . Moreover, Am 2 cannot be lower 
than ~ 3 x 10 -5 eV 2 to avoid being close to the resonance condition in the Earth and resulting 



large day /night variations of P ee . The situation is illustrated in Fig. [TJ Evidently, Nature chose 
to "tune" the mass splitting involved in solar neutrino oscillations to the density in the solar 
core! Remarkably, a completely independent reactor antineutrino measurement by KamLAND 
showed that Am 2 is indeed in this range. 
The preceding discussion assumed that mass eigenstate 3 is not involved in the evolution 
of solar neutrinos. The correction due its presence is trivially computed if we notice that the 
splitting between this state and eigenstates 1 and 2 is significantly larger than the matter 
potential even in the center of the Sun (Am 2 . tm /2E ^> \/2GpN e (0i)), so that the v e content of 
that state is always given by sin 2 #13. Repeating the arguments that led to Eq. ([T]), one gets 

Pjf > Zu = sin 4 6>i3 + cos 4 e 13 Pf e d ' 2u . (3) 

Given the bound sin 2 #13 < 0.02 from CHOOZ [12], the first term is negligibly small. The effect 
of the third state then is to multiply the two- neutrino survival probability by cos #13 . The 
resulting correction is at most 4%; this correction is basically the probability that the original 
electron neutrino "disappears" into state 3. See, e.g., [T31 [J3] for recent data analyses and 
further references. 

2. Searching for nonstandard neutrino interactions 

The impact of nonstandard neutrino interactions with matter on solar neutrino oscillations was 
discussed already in the classical paper by L. Wolfenstein [6] and subsequently elaborated on by 
many authors ( [15\ \W[ [T7] and many others) . The idea is that novel interactions due to a heavy 
vector and scalar exchange could modify the neutrino forward scattering amplitude and hence 
the oscillation Hamiltonian in matter. Regardless of their origin, at low energies relevant to 
neutrino oscillations, nonstandard interactions (NSI) are described by the effective Lagrangian 

l nsi = -2V2G F {u alp vp){e^f L YfL + e f *hYf R ) + h.c. (4) 

Here e^g (e^g) denotes the strength of the NSI between the neutrinos v of flavors a and j3 and 
the left-handed (right-handed) components of the fermions /. 

Neutrino scattering tests, like those of NuTeV [18] and CHARM [19], constrain mainly the 
NSI couplings of the muon neutrino, e.g., \e eiM \ < 10~ 3 , |e„^| < 10~ 3 — 10~ 2 . In contrast, direct 
limits on e ee , e eT , and e TT are remarkably loose, e. g., |e"" s | < 3, —0.4 < e^ R < 0.7, |e""| < 0.5, 
|e^g| < 0.5 [20]. Stronger constraints exist on the corresponding interactions involving the 
charged leptons. Those, however, are model-dependent and do not apply if the NSI come 
from the underlying operators containing the Higgs fields [21]. Here we only consider direct 
experimental bounds. 

Even with the addition of NSI the splitting Am 2 , tm /2E remains much greater than the matter 
potential anywhere along the neutrino trajectory. This means the solar neutrino analysis can 
still be reduced to two neutrino states, following the arguments of Sect. [TJ Neglecting small 
corrections of order sin #13 or higher, the corresponding matter contribution to the two-neutrino 
oscillation Hamiltonian can be written as 

H^ = ^( 1 + 611 f 12 ), where 611 = ^ " 6 TT sin 2 23 , () 
mat y/2 V £ 12 ei2 = -2e eT sm0 2 3- KJ 

The epsilons are the sums of the contributions from the matter constituents: e a /3 = 

Ylf=ude e L(3 n f/ n z- I n turn, e a a = e^jj + e3. Observe that only the vector component of 
the NSI enters the propagation effect; in contrast, the NC detection process at SNO depends 
on the axial coupling. The propagation and detection effects of the NSI are thus sensitive to 
different parameters, and the corresponding searches could be complementary. 



Eq. ([5]) shows that the flavor changing NSI effect in solar neutrino oscillations conies from 6 eT , 
while the flavor preserving NSI effect comes from both e ee and e TT . A useful parameterization is 

„ NSt ( Acos2a AT** (tin 2a \ , ^ = \e 12 \/{l + e ll} , 

<at=( 4e 2 J 0- 2a _, CO s2« h where 20 = ^(e^ (6) 

V Ae siniQ Acos2a ; A = G F n e ^[(1 + en) 2 + |e 12 | 2 ]/2 . 

The effect of a is to change the mixing angle in the medium of high density from tt/2 to ir/2 — a. 
The angle <j) (related to the phase of e er ) is a source of CP violation. Solar neutrino experiments, 
just like terrestrial beam experiments |22[ 123]. are sensitive to its effects |24j . while atmospheric 
neutrinos are not J25[ [26] . 
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Figure 2. The electron neutrino survival probability (left) and day/night asymmetry (right) 
for Am 2 = 7 x 10 -5 eV 2 , tan 2 9 = 0.4 and several representative values of the NSI parameters: 
(1) eft = e d n = e\ 2 = ef 2 = 0; (2) eft = e d u = -0.008, e u l2 = e{ 2 = -0.06; (3) eft = e d n = -0.044, 
e « 2 = e f 2 = 0.14; (4) = e{ x = -0.044, e^ 2 = ef 2 = -0.14. Reproduced from [S]- 



The main effects of NSI on P ee are as follows [24 : (i) the low-energy limit stays the 
same (vacuum oscillations); (ii) the high-energy limit changes, according to Eq. ([1]), P ee — > 
sin 2 a cos 2 6 + cos 2 a sin 2 0; (iii) at intermediate energies, the transition from vacuum to matter 
dominated regime can shift in energy, with changing A, and can become more or less abrupt, 
with changing a and <j). The nonadiabatic regime occurs when 9 — ► a, rather than 9^0. Also, 
very importantly, the day/night effect can change with all three parameters. In particular, it 
becomes small either as A — > [27, 28J or as a — > 9 [24J. Thus, the LMA-0 region that is normally 
excluded by the non-observation of day/night asymmetry may become allowed |24 t 127 1 [28]. 

Fig. [2] illustrates the impact of the NSI on P ee and the day/night asymmetry. Curve 3 gives 
an example of parameters that can already be excluded by the current data. Curve 4 illustrates 
the suppression of the Earth effect described above. For technical details, including approximate 
analytical expressions for P ee and day/night asymmetry, see |24| . 

The solar neutrino analysis of NSI cannot be done in isolation: the same NSI can also be 
probed with atmospheric neutrinos. Indeed, on general grounds, one expects the atmospheric 
neutrinos - particularly the high energy ones for which nonstandard matter effects can dominate 
over the vacuum oscillation effects - to be a very sensitive probe of NSI. Early two-neutrino 
(up, v T ) numerical studies [29] yielded e^ r < 0.08 — 0.12 and e TT < 0.2 0. Clearly, these are very 
strong bounds; if they were to extend to e eT , the large NSI effects on solar neutrinos discussed 
above would be excluded. It turns out, however, that this is not the case: when the analysis is 



Notice the difference in normalization: our e's are normalized per electron, while |29| gives e's per d quark. 




Figure 3. Left panel: A 2-D section (e ee = —0.15) of the allowed region of the NSI parameters 
(shaded). We assumed A?ti,q = and #13 = 0, and marginalized over 6 and Am 2 . The dashed 
contours indicate our analytical predictions. See text for details. Right panel: The effect of the 
NSI on the allowed region and best-fit values of the oscillation parameters. From [25J. 

properly extended to three flavors, one finds that very large values of both e eT and e TT are still 
allowed by the data [25]. This is illustrated in Fig. [3] (left panel), which shows that NSI with 
strengths comparable to the Standard Model interactions can be compatible with all atmospheric 
data. It must be noted that the compatibility is achieved as a result of adjusting the vacuum 
oscillation parameters: large NSI imply a smaller mixing angle and larger Am^ im , as can be see 
in the right panel of Fig. [3l 

The addition of the K2K data helps constrain the allowed NSI region somewhat [26]. While 
the addition of the first data from MINOS brings no further improvement [30j, the future high- 
statistics MINOS dataset will be a very valuable probe of this parameter space [30| . 

3. Searching for neutrino transition moments 

The idea that solar neutrinos could be affected by the neutrino spin precession (NSP) in the solar 
magnetic fields is even older [31] than the NSI idea. Remarkably, this idea - much improved 
with time [32j [331 EH ESI [361 EI] ~ remained viable for the next three decades. While, by the late 
1990's, the lack of time variations in the Super-Kamiokande data gave strong evidence against 
large NSP in the solar convective zone, NSP in the radiative zone continued to give a good fit 
to all solar data [38j . 

Even after the confirmation of the LMA oscillation solution by KamLAND [IT], the possibility 
of the NSP happening at a subdominant level remains of great interest, as a probe of the neutrino 
electromagnetic properties and, at the same time, of the magnetic fields in the solar interior. 
NSP coupled with flavor oscillations could lead to conversion v e — > v e in the Sun, on which 
recently, KamLAND [39J reported an upper bound. It is very important to understand what 
this bound implies for the neutrino magnetic (transition) moment and how it compares with 
other available bounds on the neutrino transition moment. 

The laboratory bounds on the neutrino magnetic (transition) moment come from measuring 
the cross sections of ve~ or ve~ scattering in nearly forward direction. The recent bound for the 
interaction involving the Majorana electron antineutrino is 2/j, e p < 0.9 x 10~ 10 /i£ at the 90% 
confidence level |43j . where = e/ (2m e ) is the Bohr magneton (m e is the electron mass, e is its 
charge). Stronger bounds, \x < 3x 10~ 12 /iB, exist from astrophysical considerations, particularly 
from the study of red giant populations in globular clusters [44]. Larger values of the transition 
moment would provide an additional cooling mechanism - via plasmon decay to vu - for the 
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Figure 4. The regions of the oscillation parameter space where one may expect the electron 
antineutrino flux above the KamLAND bound 39] (three different shadings correspond to three 
different normalizations of the magnetic field, up to the upper bound [40J). An optimistic value 
of the transition moment, /i = 1 x 10 _11 //b, was taken. For comparison, the region allowed by 
the combined analysis of the KamLAND and solar neutrino data [H] is also shown. From |42j . 



red giant core and change the core mass at helium flash beyond what is observationally allowed. 

An important theoretical consideration is that the magnetic moment operators will radiatively 
generate the neutrino mass. Requiring that the corresponding contribution to the neutrino 
masses be not much greater than their observed values, one obtains a model- independent 
"naturalness" bound on fi a b. For Dirac neutrinos, the bound obtained in this way turns out 
to be very stringent, fj, a t, < W~ 14 /j,b [IS]. However, very importantly, for Majorana neutrinos 
the bound is much weaker, only fi a b < lCT 10 /i£ (MJHTJ, owing to the different flavor symmetry 
properties of the mass and the transition moment operators [IE]. In fact, explicit models 
exploiting these very symmetry properties were discussed many years ago (e.g., [Ml [50] El). 

It turns out that for the measured LMA oscillation parameters NSP in the radiative zone 
cannot produce the v e flux above the KamLAND bound. This is illustrated in Fig. [H This is a 
remarkable example that knowing neutrino oscillation parameters precisely can be very valuable: 
the answer would qualitatively change if the mixing angle were 20° instead of 30°. 

For NSP in the convective zone, the analysis is very different, though in the end the conclusion 
is similar: one should not have expected the flux of v e in excess of the published KamLAND 
bound. Put another way, the bound on the neutrino transition moment from the KamLAND 
bound is comparable to the direct laboratory and "naturalness" bounds, but still weaker than 
that from analysis of the red giant cooling. An updated analysis of a larger KamLAND dataset 
is needed. The reader is referred to [12] for details and further references. 



2 I thank R. N. Mohapatra for bringing the last two references to my attention. 
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